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Abstract 

In our previous papers we were interested in making a reconstruction 
of quantum mechanics according to classical mechanics. In this paper we 
suspend this program for a while and turn our attention to a theme in the 
frontier of quantum mechanics itself — that is, the formation of operators. 
We then investigate all the subtleties involved in forming operators from 
their classical counterparts. We show, using the formalism of quantum 
phase-space distributions, that our formation method, which is equivalent 
to Weyl's rule, gives the correct answer. Since this method implies that 
eigenstates are not dispersion-free we argue for modifications in the ortho- 
dox view. Many properties of the quantum phase-space distributions are 
also investigated and discussed in the realm of our classical approach. We 
then strengthen the conclusions of our previous papers that quantum me- 
chanics is merely an extremely good approximation of classical statistical 
mechanics performed upon the configuration space. 

1 The Formation of Quantum Operators 

The formation problem of quantum mechanical operators has already been 
treated by a number of authors |], |[ g, g| . The problem resides basically 
in forming, from a given classical function f(q,p) of the commuting generalized 
coordinates and momenta, the related quantum mechanical operator 0[f(p, q)]. 

This problem, as expressed in ShewelPs review article ||, either cannot be 
solved in a non-ambiguous way (e.g. von Neumann's and Dirac's rules) or it 
does not give the result expected by the orthodox epistemology of quantum 
mechanics, viz. Weyl's and Revier's rules. 
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As an example, the operators related to the classical function f(q,p) = q 2 p 

are 

0[q 2 p 2 } = 0[q] 2 0[p} 2 - 2iTlO[q]0[p] - U 2 /4 



von Neumann: j Q[M = 0[g?0[p? _ 2iW[q]0 \p] - h 2 (1) 
and 

Dirac- / 0[q 2 P 2 ]=0[q?0\p] 2 -2ihO[q]0\p]-h 2 /Z 

\ 0[q 2 p 2 } = 0[q] 2 0[p] 2 - 2ihO[q]0\p] - 2% 2 /3 1 ' 

which are ambiguous and so, contrary to the postulates of quantum mechanics — 
each 'observable' must be assigned to one and only one operator. 
For Weyl's rule we have 

0[q 2 p 2 } = 0[q] 2 0[p] 2 - 2ihO[q]0[p] - k 2 /2 (3) 

which is unambiguous but has the drawback that, for the harmonic oscillator 
problem with the hamiltonian given by 

H= l -{p 2 +q 2 ), (4) 

this rule gives (figure I) 

0[H 2 ] ={0[H]} 2 + h 2 /A (5) 

which is equivalent to saying that it predicts energy dispersions for the energy 
eigenstates (that is also contrary to the orthodox epistemological considerations 
of quantum mechanics). In this example, equation (|^) says that the dispersion 
will be given by 

(AE) 2 =< 0[H 2 }- < 0[H] » 2 = h 2 /4, (6) 

for all energy levels. 

Weyl's rule has also the disadvantage that some functions of classical con- 
stants of motion are transformed into operators which do not commute with the 
hamiltonian and are not constants of motion in the quantum mechanical side. 
This may be shown using the hamiltonian || 

H=p 2 /2 + q y4, (7) 

for which, Weyl's rule gives 

0[H 2 } = {0[H}} 2 + 3h 2 q 2 /A, (8) 

implying that 0[H 2 ] does not represent a constant of motion, since it does not 
commute with H . 

In short, the acceptance of Weyl's rule implies that we do not accept the 
first von Neumann's rule || 

if 0[A] = A => 0[f(A)} = /(A). (9) 
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It is possible to show that the symmetrization rules may be reduced to Weyl's 
rule and present the same problems. 

It is the aim of the next section to show, using our formalism^ |, |, [HJ [ll], 
[l2], |l4j why ambiguities arrive and to throw some confidence upon Weyl's 
rule, which is equivalent to our own. 

In the third section, we show that this result is in accordance with the results 
obtained using the classical phase-space functions derived using our previously 
defined formalism @, |, §, [h], [l|, g || . 



Some of the main characteristics of these functions will be dealt with in the 
fourth section by means of some illustrative examples that will throw some light 
upon the discussion. 



2 The Ambiguity Problem 

It has been already shown that the ambiguity related to von Neumann's and 
Dirac's rules (to cite but a few) are related with the attempt of mapping a 
commutative ring into a non-commutative one g. We will show this using our 
analytical method. 

From our previous work one might see that the Schrodinger equation of 
quantum mechanics is derived from the classical Liouville equation using the 
Infinitesimal Wigner-Moyal Transformation 



Q- Y' q + Y' t ) = J F (^ t ) eip5qd P' 



(10) 



where F(q,p;t) is the classical phase-space probability density function and 
p (q, Sq; t) is the quantum density function (also called density matrix) and where 
Sq is an infinitesimal increment taken over the coordinate q-axis. 

To derive Schrondiger equation it is necessary to impose upon the density 
function the functional form 

,(,-|,, + | ;t )= V ,t( 5 -|;^( g+ § ;i ) (11) 

and to expand this product with 

iP(q;t) = R(q;t)e lS ^/ h (12) 

in powers of the infinitesimal increment Sq. We then take the density ([ll]) into 
its equation and, since this increment is infinitesimal, we keep only the zeroth 
and first order resulting terms and show that the former implies the continuity 
equation while the later implies Schrodinger equation. 

It is noteworthy that in this procedure we have lost one of the independent 
variables, viz. the infinitesimal Sq related to p. This makes our problem much 
simpler but it has the disadvantage of introducing the non-commutativity of 
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the canonically conjugated quantities. The explanation why this formalism 
is rather triumphant is related with the fact that the error introduced with 
non-commutativity is of the order of Planck's constant and so, the method is 
expected to be efficient to within this approximation (which is rather good). 

The construction of operators is a direct result of the definition ([lO]). Since 
the mean values of a classical function f(q,p) are given by 



< f{q,p) > 



f(q,p)F(q,p] t)dqdp, 



(13) 



it is automatic that the mean of the operator related with this function is given 

by@ 

p(q-^-,q + %t)dx (u) 



< 0[f(q,p)] >= / lim / (q,-ih 

J Sq~*0 \ 
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d(Sq) 



giving the same values of expression (p^). 



The funct 



ion / (q 



-ih 



is unambiguous defined since q and 5q are 



independent variables. 



Now, giving the amplitudes ip (q; t), we can form p [q — 4?, <Z + ^2 it) accord 



So. 



ing to expression ( |11| ) and expand these amplitudes in a Taylor series around q 
to get, until second order in Sq 



< 0[f(q,p)] >= / lim / q,-ih 



1 

'2! 



Sq^O 
2 



d 



d(5q) 



Sq 



I dtp dtp^ 



dq 2 dq 2 dq dq 



dq 
o(Sq 3 ) 



dq 



(15) 



and then apply to each term of the expansion the operator / defined in ( |14[) 
and take the limit as indicated. 

The procedure defined above is unambiguous and, it might be easily shown, 
gives the same results of Weyl's rule. We have, for example, in the case of 
f(q,p) — q 2 p 2 , the operator 

d 2 



f q, -ih 



d 



= lim q —fi 

Sq^O 1 



d(Sqf 



(16) 



giving, in terms of the amplitudes, the operator 

f = q 2 f-2ihqp-h 2 /2. (17) 
However, in expression ([li]), we could have considered the term q 2 as part 
%,q+ S -§ 



are calculating 
<0[f(q,p)] >-- 



of the function p [q — ^-,q + ^;t); in the example above this means that we 



lim 

Sq^O 



, d 2 

d{5qf 



q p\q 



Sq Sq 



dq, (18) 
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which is allowed since q and Sq are independent variables. In this case we will 
expand the function 



Sq Sq 



q p l q - y,<?- 



t) =<?V [q-~-^-,t)ip[q + 



Sq 



Sq 



(19) 



and there are three possibilities of doing that 



W (q - s -§ ; t) }i> (9 + s 4 ; «) ? 2 p 2 - 2z^5p - ?* 2 
V^t {<z 2 V (<z —> ? 2 £ 2 - 2*n$p - ft 2 

{#t (g - §£; i) (? - 4?; *) } — > q 2 p 2 - ^kqp 

or we could make the expansion of 
d 2 



(20) 



d{Sqf 



} 



(21) 



giving two other possibilities. 

All these results are different from the previous one (|l7| ) because the in- 
finitesimal increment is taken over the coordinate g-axis and taking the limit 
Sq — > implies making one variable to disappear. 

Another way to see that is making the coordinate transformation 



, &q , 
y = q + Y>y 



q 7t 



6q 
2 



giving 



d 



d(Sq) \dy dy 



We then have 

00 ^ 

p(y',y,t) = y^ —, 



dy' 



(22) 
(23) 

p(y' ,y ;t) (24) 



where it is implicit that we first perform the derivatives and then compute their 
values at 

yo = y'o = q (25) 

Using expression ( p3f ) we get for the mean value of the classical phase-space 
monomial represented by q N p M the corresponding expression 



< q N v M >={-ih) M J { l±p- 

00 1 











71 = 



d d 



dy dy' 



p(yo,yo;t)d 



M 



y + y' 



(26) 
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Now, since y' and y are independent coordinates we might change the order 
of some of the terms in the above expression which are of the form 

■ (27 » 

but when we calculate the limit we see that, in this limit, 

y = y' = q (28) 

and are not independent quantities anymore (they are, indeed, equal). This 
limiting process is exactly the one we have used, as described above, to derive 
the Schrodinger equation and is also, as already shownjij the one responsible 
for the introduction of non-commutativity of the dynamic operators. 

We are not able to choose, at this place, one or another of the procedures as 
shown in the examples (|l^,^) and (|2l| ) but the claims for the procedure leading 
to Weyl's rule are obvious. This is because it is not clear why we will expand 
that part of the operator (we are trying to find an expression for) as part of the 
density function; indeed, in expression (|l6|), for example, the term q 2 refers to 
an operator while the density function is indexed by q, a coordinate variable. 

We must be aware however, of the consequences of this choice. One of them, 
for example, will be to predict, for the harmonic oscillator example, a dispersion 
in the energy according to expression (j^) of 

A£=|, (29) 

a result to which we will soon return. 

In the next section, we will show that these claims for Weyl's rule is strength- 
ened when we return, by means of expression (|To|), to the classical phase-space 
formalism. We will also discuss some of the properties of these functions in the 
realm of the present developments. 



3 Classical Phase-Space Distributions 

Distributions of probability defined upon classical phase-space were first pro- 
posed by Wigner |l5|| in connection with thermodynamics considerations. Lat- 
ter, the theme was further developed by Moyal|l6|| and GroenewoldQ. In the 
late sixties, the problem began again to deserve some attention mainly in con- 
nection with problems on the coherent properties of the electromagnetic field 
on the theory of laser and on non-linear processes jl], [|. 

More recently, we have seen the revival of the problem of writing quantum 
mechanics using these distributions and the notion of stochastic phase-space 

f|[§ §§ n n, n §, n n u. 

All these latter developments are rather technical and pay more attention to 
formal aspects of these distributions than to the epistemological considerations 
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that could follow from these developments. Such considerations were indeed 
well formulated in Moyal's paper where the epistemological deficiencies related 
to the tentative of linking classical and quantum mechanics by writing the later 
in classical phase-space were pointed out. 

The major problem with the distributions F(q,p; t), derived from expressions 
nearly similar to (|h]) and defined over classical phase-space — as recognized by 
all the cited authors — , is the fact that these functions might take negative 
values; they are then called quasi-probability functions. This fact is usually in- 
terpreted as meaning our necessarily introduction of an error when passing from 
a fundamental theory (quantum mechanics with non-commuting 'observables') 
to an approximate theory (classical statistical mechanics with commuting dy- 
namic variables). 

In the context of the present series of papers this conclusion has to be 
reversed. We have already shown that it is possible to derive quantum me- 
chanics from classical mechanics without introducing any non-classical postu- 
late. We have also shown that this process of derivation, where one of the 
phase-space coordinates (usually the momentum) disappears, introduces the 
non-commutativity when we impose that our densities have to obey equation 
(p"l|). The hierarchy of the theories is fixed by this derivation and we see that 
quantum mechanics will give good results whenever the experimental environ- 
ment is not able to uncover Nature with techniques where the dispersion product 
is given by 

AqAp < |, (30) 

where q and p are the generalized canonically conjugate coordinate and momen- 
tum. That is to say: the dispersion relations do not impose any constraint upon 
the behavior of Nature but only upon our capacity of describing Nature by means 
of quantum theory. 

Because of the smallness of Planck's constant, we might conclude that this 
technique is rather fruitful and the error introduced by the quantum formalism 
is far from being detectable with nowadays laboratory technology. 

Then, as we have already shown, the passage from a description based upon 
the classical density function to a description based upon the density 'matrix' 
p(<7, Sq) is exact but we still have two independent variables and, as in the 
classical phase-space, no dispersions 0. When passing from densities p{q,5q) to 
a description using the amplitudes ip(q] t), the dispersion relation becomes 

AgAp>|, (31) 

but the problem becomes rather tractable and, as the error introduced is ab- 
surdly small, we might conclude that more was gained than was lost. 

Now, when trying to return to classical phase-space — not by inverting the 
Fourier Transform ( |To| ) but by finding the function F(q,p;t) that, when inte- 
grated will result in the desired p(q, <5<z)f| — we shall admit that the introduced 
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error related with the dispersion relation and with the imposition ( p"i"| ) will prop- 
agate upward and will necessarily be incorporated into our densities F(q,p;t). 
This fact might be easily verified from the usual expression for the classical 
density 

F(q,p;t) = e -^ hd2 / dpdq W {q)4>{py m,h \ , (32) 

where ip(q) and 4>{p) are the position and momentum eigenfunctions. This 
expression does indeed incorporate Heisemberg's inequality. 

The fact that such functions are not positive definite might be taken as 
an indication that the process of quantization — equation ([IT]) and the limit 
Sx — > — indeed introduces an error. 

Another problem with these distributions is the fact that, although they give 
the correct values for the energies mean values and the dispersion relations for q 
and p, they also predict a dispersion in the energy eigenstates — something rather 
contrary to the orthodox interpretation where eigenstates are dispersionless. 

This is the point we wanted to reach that will help us with the discussions 
about which method of operator construction to choose. 

In the next section we will develop some examples of classical phase-space 
distributions for quantum systems and will study more deeply the harmonic 
oscillator problem to stress some relevant points for our latter discussion. 



4 Some Examples: The Harmonic Oscillator 

We begin by giving some obvious examples: 
(i) Free Particle: 

i/>(x,t) = e lkx p(x, Sx) = e lkSx =► F(x,p;t) = 5 D (k - p/h) (33) 

where Sd(x) is the Dirac's delta function. This density p(x, 8x) shall be com- 
pared to the one obtained when passing from the classical phase-space to the 
'quantum' (x, Sx) configuration space. Indeed, we have already shownQ that, 
for a dispersion free ensemble of systems with one free particle, where the clas- 
sical density function is given by 

F(x,p;t)=S D [x-xo(t)]S D \p-p {t)], (34) 

the correct density function is given by 

p(x,Sx) = S D [x - x (t)]e ik5x (35) 

and represents an ensemble of systems with one localized free particle. Since 
this density function represents a dispersion free ensemble it cannot be expressed 
within quantum theory based upon Schrodinger's equation. We thus clearly see 
that, when passing from the quantum to the classical, this localization cannot 
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be retrieved and we are left with expression ( J33| ) which has no information about 
the position of the particle — this information was lost, 
(ii) Minimum Free Wave-Packet: 



ip(x,t) 



„ — x /2a -\-ik x 



p(x,Sx) = e -^ 2 e -(Sxf /io?+zk Sx ^ 



where a — Ax\/2. We then have 

F(x,p;t) = Ne -F/2 a 2 +(P^k fa 2 /n 2 }^ 

which implies that Ap = hsplja and so 

AxAp = h/2. 

We now want to develop another example more deeply. 

(iii)The Harmonic Oscillator: The probability amplitudes are given by 



(36) 
(37) 
(38) 



1 



2"n!7r 1 /2 



1/2 



H n (x)e 2 X \ 



(39) 



where H n (x) is the n-th Hermite polynomial. These amplitudes are solution of 
the Schrodinger equation 

Ei> n {x) = g [f + Q 2 ] ipn{x) = E n ij} n {x) = (n + l/2)hi/j n (x), (40) 

where E is the hamiltonian operator and E n its eigenvalues. 

Using expression ( |Io| ) we might explicitly calculate F(x,p;t) to get the nor- 
malized functions 



n = 
n = 1 

n = 2 
n = 3 

where 



F (x,p;t) 
Fi(x,p;t) 

F2(x,p;t) 

Fz{x,p;t) 



-H 



-H 



2 
37rfi 



(H-l) 



2(H- 2) 3 + (H — 2) 2 - 3(i? - 2) 



2 /t2 



-if 



(41) 
(42) 



In figures II, III and IV we plot the choices n = 0, 3 and 10 respectively. It 
is clear that, apart from the ground state, we have the appearing of negative 
probability densities. 

As amply noticed in the literature, we get non-negative densities when we 
integrate in one or other phase-space coordinates giving the lateral probabilities 



F^x; t),F 2 (p; t) > Vx and Vp 



(43) 
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which are all in accordance with the quantum mechanical results. 

With the density functions F(x,p;t) at hand we might calculate the mean 
energies, the Heisembcrg dispersion relations and the mean quadratic deviation 
(the dispersion) of these energies. We calculated these values for the first n = 
0, .., 7 functions and quoted the results in table I. 

Looking at this table we can see that the energy mean values and the Heisem- 
berg dispersion relations are given by 

< E n >= (n + l/2)h and AxAp = (n + l/2)h (44) 

as expected Jl7|] . For the energy mean quadratic deviation we get 

AE = \ (45) 

for all states. 

We now return to expression (||) and see that result ([45| ) agrees with the 
way we define operators using our formalism (and also with Weyl's rules), ft is 
important to stress that, when finding F(x,p;t) from the quantum mechanical 
amplitudes we have not made any approximations — we have just performed the 
'inversion' indicated by (|Io|). The value of the energy dispersion shall agree with 
the quantum mechanical one as in the cases of the mean energy and Heisemberg's 
relations (||). 

We then must abandon the first von Neumann's rule|| represented in ex- 
pression (||) and also conclude that even eigenstates of the Schrodinger equation 
are not dispersion-free. 

The results above are good indications that our operator construction proce- 
dure, and also Weyl's rule, is the correct one to be used in the formal structure 
of quantum mechanics. 

Another consequence of expression (pOf), which may be written generally as 

AE = ^uj, (46) 

where lo is the frequency of the harmonic oscilator (we have made equal to 
unity in the previous development), is the conclusion, based on the Heisemberg's 
energy-time dispersion relation, that the dispersion in time must be given by 

1 T 

At=- = -, (47) 
lo 2ir 

where T is the period of oscilation. 

We also note that, considering expressions ( [44] ) and (f|6|), the correct energies 
(within the mean quadratic deviations) of the harmonic oscilator problem shall 
be written as 

E„ =< E„ > ±AE = {n +-)h ± - , (48) 

which throws some light upon the meaning of the vacuum energy fluctuation 
(we will return to this problem in a latter work). 
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5 Constants of Motion 



The acceptance of our rules of operator formation leads to the problem eluded 
in (^) where functions of constants of motion are not necessarily constants of 
motion — since, now, some function of the hamiltonian operator, say f{H), might 
not even commute with H. 

We might see this in a straightforward way if we look for the equations 
satisfied by the density functions F(q,p;t). If we write F®(q,p;t) for the func- 
tion obtained from the quantum mechanical amplitudes, as above, and write 
F c (q,p;t) for the classical distribution, we have, by hypothesis 

— + {H,F c } = (49) 

while, it is possible to show]Tr|, the equation obeyed by function F®(q,p;t) is 
given by 

f)F® 9 DBF) 

-ar + r in 2^'^^^^ ; *)=°- (50) 

This makes clear our point of view. We began with equation ( (49| ) and end 
with another different equation (^). This is another argument for concluding 
that we shall not have 

F c (q,p;t) = FQ(q,p;t). (51) 

In some special cases (the harmonic oscillator is an example) equation (^) is 
equivalent to equation (^0|) . Even in such cases we shall not have (^lj) since we 
have postulated, in passing from the former to the later, one specific functional 
form for our densities and also a dispersion relation given by (^l|). 

It is also important to stress that the 'classical limit' cannot be generally 
attained by taking the limit h — > as usually stated in the literature. This was 
also shown elsewhere jl8| . 

The important result ( |50j ) means exactly what we have said above. Indeed, 
even if some quantity does commute with the Hamiltonian operator, some func- 
tional of it (e.g. its powers) will not generally satisfy this requirement and 
will not be a constant of motion |[(|. This is another result that indicates that 
expression ( |50| ) must be only an approximation. This was already exemplified 
above in expressions (Q) and (@), where we have noted that the square of the 
hamiltonian does not commute with the hamiltonian itself and cannot be a 
constant of motion. 

This same result might be exemplified from Dirac's rule. Since this rule uses 
precisely the substitution of the classical Poisson brackets by the commutator || 
and does not produce, in general, unambiguous operators, the one-to-one cor- 
respondence of classical and quantum-mechanical constants of motion generally 
fails. 
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6 Conclusions 



Until now we have been concerned with reconstructing quantum mechanics from 
classical mechanics || § g |l0[ 0, [l§ [H| 0. This task was taken from the 
mathematical and epistemological points of view and there are, still, much work 
awaiting — for example, the reconstruction of quantum field theory. 

This paper, however, has undertaken a different approach. The theme we 
have chosen to approach here lies in the frontier of applicability of the quantum 
formalism itself. 

Indeed, one of the most fundamental axioms of the orthodox interpretation 
of quantum theory is the one-to-one correspondence between classical functions 
(to be quantized) and quantum operators (observables) . We have thus avoided 
to contradict this postulate since it comes naturally from our own derivation of 
quantum mechanics. 

We were then led to adopt Weyl's procedure, or the one we have derived from 
our own theory, as the correct one for the formation of quantum mechanical 
operators. The most striking consequence of our choice was to be left with 
energy eigenstates presenting energy dispersions. 

We then have shown that the appearance of these dispersions are pertinent 
to the formalism since we derived them from the quantum phase-space distri- 
butions F^(x,p;t), defined above, which keep the same physical content (they 
are only Fourier transforms of the quantum density function) of the quantum 
distributions. 

The unavoidable appearance of energy eigenstates with energy dispersion 
points in the direction of an interpretation of these states very different from 
the one given by the orthodox epistemological approach. We must abandon, 
for example, the idea that when performing a measurement upon an energy 
eigenstate we do not disturb this energy value and the next measurement will 
observe the same value for this variable. 

Concerning the formalism, the adoption of our procedure of operator forma- 
tion leads to the physically unacceptable result according to which functions of 
constants of motion are not themselves constants of motion. It is interesting to 
note, however, that this pathological behavior comes from the attempt to map 
a commutative ring into a non-commutative one, and is of the order of, at least, 
the square of Planck's constant; since the magnitude of this error is highly be- 
yond experimental precision, we may still expect this pathology not to be easily 
empirically detectable — this, of course, does not preclude us from studying the 
theme, since we are interested here in placing the quantum theory on the most 
solid ground. All these results might be also derived from the quasi-Liouville 
equation (|5(]) obeyed by the quantum phase-space distributions above. 
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n 


energy 


AqAp 


AE 





1/2 


1/2 


1/2 


1 


3/2 


3/2 


1/2 


2 


5/2 


5/2 


1/2 


3 


7/2 


7/2 


1/2 


4 


9/2 


9/2 


1/2 


5 


11/2 


11/2 


1/2 


6 


13/2 


13/2 


1/2 


7 


15/2 


15/2 


1/2 



Tabic 1: Mean energies and dispersions for the harmonic oscilator derived from 
the quantum phase-space distributions. 
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Figure 1: Graph showing the fail of von Neumann's process of quantization. 
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